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Abstract— We report that an exact analytical solution 
for the modal electromagnetic fields in a conducting 
rectangular waveguide filled homogeneously with an un- 
rotated biaxially anisotropic magneto-dielectric material 
can be obtained directly by separation of variables, 
without recourse to numerical methods. It is found that 
higher order modes are generally hybrid non-TE/TM type 
modes. It is shown for strongly anisotropic media that 
there exist pseudo-evanescent complex conjugate mode 
pairs. This is the first time complex conjugate mode pairs 
have been directly demonstrated analytically in waveguide 
completely filled with magneto-dielectric material. Lastly, 
it is demonstrated that no alternative mode sets exist. The 
analytical solution given in this paper has not yet been 
reported in the literature, and could prove useful in 
guiding efforts to develop  Nicolson-Ross type 
measurements of non-ideal anisotropic behavior exhibited 
by passive, micro-structured, magneto-dielectric 
materials. The solution also adds to the portfolio of exact 
analytical solutions in applied electromagnetic theory, 
which can be utilized to test the veracity of general 
purpose numerical codes. 


Index Terms— Complex modes, anisotropic materials, 
magnetic materials, electromagnetic theory. 


I. INTRODUCTION 


In this report, an exact analytic solution is given for the 
electromagnetic modes supported by a conducting rectangular 
waveguide filled homogeneously with a composite magneto- 
dielectric material described by permittivity and permeability 
tensors. Such materials can be realized in several ways: by 
structuring of multiple different material phases [1-3]; by 
treatment of soft bulk ceramic polycrystalline ferrites in strong 
magnetic fields during crystallization [4]; by formation of 
bulk ferrites in the form of a single large crystal [4]; or by 
constructing thin-film composites with preferentially aligned 
single-domain ferromagnetic metallic alloy flakes, strips, or 
rods [5-6]. In developing such composite materials, the crystal 
structure or orientation of the inclusions is controlled to 
enhance the magnetic or electric properties of the material, but 
often with the side effect of making the material’s properties 
directionally dependent, or anisotropic. 


The mathematical boundary value problem for modeling 
the magneto-dielectric-waveguide system was considered by 
Meng et al [7]. We note that in [7] an analytical solution is 
presented. However, their formulas are only valid for the TEon 
and TEmo modes. The formulation presented in [7] yields field 
solutions, that for all modes beside the TEon and TEmo modes, 
do not satisfy Maxwell’s equations. Hence, their formulation 
does not provide accurate expressions for the remaining 
higher order modes. The subsequent numerical analysis in [7] 
was only applied in the special case of uniaxial anisotropy, 
and in this case, their field solutions reduce to forms that do 
correctly satisfy Maxwell’s equations. 


Here, we give a derivation of the complete mode set for the 
case of biaxial anisotropy, including the remaining higher 
order modes. Such theoretical information is potentially of 
importance for developing Nicolson-Ross type measurement 
procedures for the empirical investigation of the non-ideal 
anisotropic characteristics of microwave magnetic materials 
[8]. 

The mathematical procedure is essentially separation of 
variables. However, the anisotropy of the magneto-dielectric 
filling necessitates the solution of an eigenvalue-eigenvector 
problem for the propagation constants and field components 
of the individual modes, which adds considerably to the 
complexity and richness of the observant electromagnetic 
phenomena. 





Fig. 1. Geometry of the waveguide. 


II. THE RECTANGULAR WAVEGUIDE FILLED WITH BIAXIAL 
MAGNETO-DIELECTRIC: SEPARATION OF VARIABLES 
SOLUTION 


We presume that the walls of the rectangular waveguide are 
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perfectly conducting. The geometry of the waveguide is ðH, j 0? E; 3’ Ey 
illustrated in Fig.1. The waveguide is filled homogeneously Ao uzo UzNoko Exlgko Ex + >> ay? ~ axdy » a4 


with a biaxially anisotropic magneto-dielectric material 


aligned with the Cartesian coordinate axes such that it can be where no = ./Uy/Eq and ky = w,/UpE9 . We define the state 


described by the constitutive tensors vector for the electromagnetic system to be the column vector 


= 0 0 given by: 
x 
Z= i e 0] , (1) wp=[Ex Ey Hy Hy)". (15) 
Ez 
Inspection of (11)-(14) enables us to reformulate Maxwell’s 
u 0 0 equations as the linear operator equation, given in state-space 
a=|0 m 0], (2) representation as: 
0 0 KL 


W iig, (16) 


where L is the linear matrix operator 


Which, assuming e/“* time dependence, enter into Maxwell’s 
curl equations as follows: 


VX E=—jopoh: A; 





















































(3) 
Sats 0 0 Ly3 Ly, 
VX H=jwéé'E. (4) i 0 DS sae Li i 
P| Lz Lz 0 0?’ CD 
It is convenient to represent the tensor equations (3) and (4) In, Laz 0 0 
explicitly in terms of their scalar components in Cartesian . 
coordinates with elements 
3 2 1 
ðE, Ey , TAE ie (18) 
ay a VEES Si 13 ezko 0xðy ' 
dE, OE. jno a” 
x B L = ky +=), 19 
az ago —JWUobyHy , (6) 47 zk, (ie o + 7) (19) 
ðE, OE. j a? 
y X . Jno 2 
—- = —jwHokHzH; , (7) z> (x e,k +5) ‘ (20) 
ax dy ZAAZ 23 Ezko x“z™*0 ay2 
= E, 8 _ Jo 
ay az 5 jwEo£xEx (8) L24 > Rory (21) 
OH, OH, . j 0? 
az Ox Ry ls, 9 -— a, (22) 
MzNoko Axdy 
OH, OH. 
y xs ; 2 
nAn = jWEgEzE; . (10) j 2 0 
0 0 L —— k =], 23 
x y 32 > usnoko EyHzKo Tayo (23) 
In a source free region, all six field components are not i 22 
independent, and any two may be eliminated. We will hes ( Exhzko? +5) (24) 
proceed by solving (7) and (10) for H,,E,, and then HzNoko dy 
substituting into the remaining equations to obtain four i 22 
coupled partial differential equations involving the transverse Lo > AME ind, f (25) 
field components, MzNoko Axdy 
We presume variables can be separated, and suppose that the 
OE a? Hy H a? H modal fields take the form 
4 -Ezko Ax? Ax y E, = E, cos(k,x) sin(kyy) 7 (26) 
ðE, _ jno > 07H, 07H, =~. 
y ky? Hy = 12 Ey = Ey, sin(k,.x) cos(k ‘ (27) 
i ele + D Iy (12) y T Ëy x (kyy) 
H, = H,sin(k,x) cos(k > 28 
dz uznoko |2 0” = ~ Axdy Hy = Hy cos(kyx) sin(kyy) , (29) 
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where, for any positive whole numbers m and n, k, = mm/a 
and ky = ni/b. The particular choices of cosine or sine in 
these expressions are made to satisfy the boundary conditions 
at the conducting walls of the waveguide. We refer to 
Ey, Ey Hees and H, as the spectral fields, which physically 
represent the weighting coefficients for the individual modes, 
and we define the spectral electromagnetic state vector to be 
the column vector 


y=- & mh Bl) (30) 


By substitution of (26-29) into (16), the spectral state vector 
can be shown to satisfy the ordinary matrix differential 
equation, 


dp 


LA UAT (31) 
dz Y 
where F is the 4x4 matrix 
= [0 Z 
r= io 4 i 32 
vO (32) 


and where we have also defined 2x2 impedance and 
admittance matrices 








— j —k,k Kye? — py Egy” 

z= a : x Ha , (33) 
Ezko [uxEzKo —ky kyky 

7 kk Eliko =k; 

Y = irl > xy ; ylz 0 x | ; (34) 
Hznoko |k," — exuzko —kyky 


and 0 as the 2x2 matrix of zeroes. Further, we define column 
vectors for the spectral transverse electric and magnetic fields: 


T=% 5]; (35) 
m=i m] (36) 


Then, (31) splits into the two separate matrix differential 
equations, 


dE, Z= 








=Z-H,, (37) 
dz t 
dH, a 
-y-E.. (38) 
dz t 
Taking the derivative of (37), we obtain 
= = 
PE See. (39) 
dz? dz 
Substitution of (38) into (39) yields 
E = 
—-7-Y-F. 40 
q ZYE. (40) 


Equation (40) is a two-equation system of second order 


differential equations. Substitution of a trial solution of the 
4v1z3 


form e*‘“7D into (40) gives 
d? — (e ep | -7-Y (e+ 777) 
dz 
tvdz (Z-Y AB) =0, 
Z‘Y-V=ADv (41) 
tViz% 


Equation (41) indicates that e 
solutions of (40) provided that ¥,/ are an eigenvector and 


are two independent 


associated eigenvalue, respectively, of the matrix ZY. Since 
Z -Y has two eigenvectors, we have four total solutions to 


(40). It follows that +a, i +a, are the four eigenvalues of 


the matrix I’. For the remainder of this paper, when we say 
eigenvalue, we unambiguously mean either A, or Az. We write 
the general solution of (40) as, 


Ge VaT + GeV uT + Ger VaT; 
+ Gev2 7, 


E, = 
(42) 


where &-ĉ; are arbitrary constants, and V7, Vz are the two 
eigenvectors of the matrix Z-Y. The transversal spectral 
magnetic field may be obtained using (37) and (42) as, 


H, = -ê me- VEZ- 7; + êh eyaz. 
Ti — & le VZ -T + GeV eZ T. 


Equations (42) and (43), taken together, form the general 


(43) 


solution to the original matrix equation (31). 


(44) 


D> 
+é3e z| a | 
—JA,Z~* + Vz 
z vz 
+é,ev222 ana 
A,LZ* + V2 


Equation (44) may then be combined with (26)-(29) to 
synthesize the field expressions for the complete mode 
spectrum supported by the waveguide filled homogeneously 
with biaxially anisotropic magneto-dielectric. The axial fields, 
if desired, may be obtained from (7) and (10). 

We note that (44) indicates two mode types, one propagates 
with propagation constant NEN , the other 4/43, and the field 
components exist in different proportions according to the 
form of the two different eigenvectors. This result is at 
variance with equations (12) and (19) of [7], which assign two 
propagation constants to a single mode type. A single mode 
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type can have only one propagation constant. This is true even 
for birefringent crystals, where birefringence results from 
coupling to two separate modes, each with different 
propagation constants, not coupling to a single mode with two 
different propagation constants. 


The modes of both types, of various orders m and n, are 
obtained from (26)-(29). Equation (44) shows that the 
waveguide supports propagation of both mode types in either 
direction, verifying the expectation that mode propagation is 
reciprocal. That is, the four solutions to (31) consist of the two 
mode types, forward and backward propagating. This is 
known not to be true for magnetic materials exhibiting a 
permanent or forced DC magnetization [9]. We could have 
also concluded immediately that propagation would be 
reciprocal since i = ff’ and = = g", which is a sufficient 
condition for the Lorentz reciprocity theorem to hold [10]. 


IJ. THE PROPAGATION CONSTANTS AND EXISTENCE OF 
COMPLEX CONJUGATE MODE PAIRS 


The solution to the eigenvalue equation (41) may be 
obtained by techniques of linear algebra. For the 2x2 
eigenvalue-eigenvector system, explicit formulas exist for the 
eigenvalues and eigenvectors. Applying these formulas leads 
to the following solutions. We obtain for the eigenvalues 

+ *) ky? = (exl + 


wiene 


2L\ez Mz 
ety ko?) + fE (6-5) + (2-2) by? + 

1/2 
2) (2-2) 


(Uxéy — Elly )ko?] Hike ky” (= £ 
where + indicates the two different eigenvalues 442 of the 
matrix Z : Y, corresponding to the two different mode types. 
A single formula cannot be given for the eigenvectors, and a 
number of special cases must be enumerated. All special cases 
can be derived by solving (41) subject to (45). 


(45) 


Case I: ky + 0, ky + 0, €yUz # EzMy: 


kyky (2 = =) 


Mae Aya — ky? - ky” + Explyko” oo 
Case II: ky # 0,ky # 0, € Mz # EzMx: 
Ana — Eka? — Sky? + Haty ko” 
V12 = ie, é T bs ) (47) 
Case II: ky # 0, ky + 0,8 x ji: 
s-a[i}eall as 


Case IV: k, = 0, ky # 0: 


4 


=. w 
Case V: k, + 0, k, = 0: 
s= [0]. o 


There is physical significance to each of these special cases. 
Cases I and II correspond to general biaxial anisotropy, and 
(46) and (47) result in the higher order hybrid modes for which 
the electric and magnetic fields in the direction of propagation 
are both nonzero. Two cases are needed to accommodate the 
possibility of partial symmetries between the permittivity and 
permeability tensors, which zero out the off-diagonal entries 
in the matrix Z : Y, changing the character of the eigenvalue 
equation (41). In the case that neither of the partial symmetries 
occur, the eigenvectors given in equations (46) and (47) are 
directly proportional, yielding the same solution. The hybrid 
modes are the modes that we are particularly interested in. We 
shall summarize the mode behavior for cases III-V before 
further investigation of these higher order hybrid modes. 


Case III is for a very special type of anisotropy where the 
permittivity and permeability matrices are directly 
proportional. In this highly contrived scenario, the matrix Z - 
Y becomes diagonal. The eigenvalues are repeating, creating 
a degeneracy requiring two eigenvectors to be assigned to a 
single eigenvalue. Equation (48) will lead to TE*/TEY mode 
sets, but other classifications could also be made by taking 
combinations of the degenerate modes. This procedure is 
discussed in the context of isotropic media in the next section. 


Cases IV and V lead to TEon and TEmo modes respectively, 
regardless of the form of the permeability or permittivity 
matrices. They are the only TE/TM modes existing in the 
general biaxially anisotropic guide. These modes can also be 
found in the paper by Meng et al, provided that we set k, = 0 
or k, = 0, and only retain the particular propagation constant 
that is found in agreement with our formula [7]. Depending 
on the dimensions of the waveguide, either the TEo: or the 
TE,0 mode is dominant in the anisotropic guide, which is no 
different than in an isotropic rectangular waveguide. For 
Cases IV and V, only Vj is given, because the fields for the 
other mode type vanish when either k, = 0 or k, = 0. 


We note that another special case occurs when biaxial 
anisotropy reduces to uniaxial anisotropy, such as when uy 
Hy and £y = £. This, however, is already contained in cases 
I-III. For uniaxial anisotropy, all modes are either TE or TM. 
The last special case is fully isotropic material, which will be 
discussed in the context of this derivation in the next section. 


Equations (45-50) reveal a complicated dependency of the 
modal fields on the permittivity and permeability tensors. 
Importantly, it appears from (45) that even if all elements of 
the permittivity and permeability matrices are real, there is, in 
general, no single well-defined cut-off frequency for each 
mode. This can be seen by setting the eigenvalues in (45) 
equal to 0, and noting that the resulting equation will not 
always have exactly one positive real root for the square of the 
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cutoff wavenumber ky’. In any situation where there is not 


exactly one real, positive root of (45) for ko’, itis not possible 
to define a single cutoff frequency in the usual way. 


To illustrate this behavior, we now consider a typical 
example. For simplicity, the material’s properties are taken as 
non-dispersive. This, of course, is not correct; especially for 
magnetic materials which exhibit various resonances 
including domain wall resonances and ferromagnetic 
resonance [4]. However, taking this shortcut will allow us to 
demonstrate the strange behavior of the propagation constant. 
We let: Ex = 2, &y = 4, €; = 15, Uy = 8, 4y = 3, Uz = 
17,a = b = .05 (meters), m = n = 1 . In Figs. 2 and 3, we let 





J = qi +jfi, and plot æi, ßi, respectively, versus 
frequency f for frequencies ranging from DC up to 2 GHz. 
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Fig. 2. Attenuation coefficient vs. frequency. 
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Fig. 3. Wavenumber vs. frequency. 
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We observe from Figs.2 and 3 that there are no well-defined 
cutoff frequencies for either mode type. The odd-looking 
kinks in the plots result from mapping the real and imaginary 
parts of a complex function of frequency onto separate plots. 
There are multiple points at which the eigenvalues are equal 
to 0, and they demarcate transition frequencies at which the 
behavior of the waveguide transitions. It can be seen that 
although a single cutoff cannot be defined, there is a minimum 
frequency at which non-attenuated propagation can occur. 
Above this frequency, the behavior of the wavenumber is 
analogous to what would occur in a waveguide filled with 
isotropic material, and the modes carry power. Below these 
frequencies, there exists an anomalous dispersion frequency 
region where complex modes propagate with attenuation, and 


exhibit backward-wave behavior. In our example, the onset of 
these complex modes is at approximately 500 MHz. It can be 
deduced from (45) that in the anomalous frequency region, the 
complex modes appear as a conjugate pair, which is necessary 
to ensure that taken as a pair the pseudo-evanescent modes 
carry no overall time-averaged real power. In our example, it 
would make little sense for only one of the m = n = 1 modes 
to exist without the other, since it would carry time-average 
power, but with attenuation, implying a loss of energy without 
a physical mechanism, as the material parameters need not, in 
theory, have an imaginary part for this behavior to manifest. 
For physically realizable fields, it follows that the presence of 
one such mode must imply the other. 


The existence and properties of such complex mode pairs 
was discussed by Rozzi et. al. for the case of lossless isotropic 
closed waveguides of arbitrary cross section, and were 
explained as another reactive mechanism of the waveguide, 
analogous to a single real mode below cutoff [11]. Indeed, in 
agreement with the requirement set forth by equation (4) of 
[11], we note from (43) and (46-47) that the transversal 
electric fields of the two modes in the anomalous frequency 
region are the complex conjugate of each other. To the 
author’s knowledge, in the context of an anisotropic 
waveguide that is closed, lossless, and reciprocal, this 
phenomenon has not been directly demonstrated before now. 


IV. COMPARISON WITH FORMULAS FOR AN ISOTROPIC 
MATERIAL: NONEXISTENCE OF ALTERNATIVE MODE SETS 


It is of interest here to show that our formulas reduce to the 
correct form when we set €y = £ = €; and Uy = Hy = Hz. 
First, in isotropic media, we note that (45) reduces to 

_ 2_4,2 2 2 
A=-k,o =k, + ky” — ke. (51) 
Equation (51) is recognized as the well-known separation, or 
dispersion, equation in isotropic media. We drop the 1,2 
subscript, because the eigenvalues (propagation constants) for 
both mode types have the same form in an isotropic 
waveguide. 


We can obtain the eigenvectors in isotropic media from 
(41), but the general character of the eigenvalue equation is 
different for repeating eigenvalues. In fact, substitution of 
(51) into (41), and simplification of the expressions for the 
matrix Z-Y yields the following equation, which places 
essentially no restriction on the eigenvectors: 


0 07 =_ 70 
lo al kz lol i 
This has, among other possible solutions, the trivial solution, 


s-a] 


(52) 


(53) 


Equation (53) shows that the eigenspace for the repeated 
eigenvalue in isotropic media is two-dimensional; spanned by 
two linearly independent basis vectors. We can rotate to other 
basis vectors that span this space and obtain alternative mode 
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sets such as those described by Harrington [12]. Rotating the 
basis vectors within this eigenspace is equivalent to taking 
combinations of the degenerate modes, that is, modes with the 
same propagation constant, to form a different set of modes. 
This can only be done when dealing with isotropic media. We 
note that in the biaxially anisotropic magneto-dielectric 
waveguide, excluding case III described above, eigenvalues 
do not repeat, the propagation constants are unique, and 
therefore alternative mode sets do not exist. 


The usual TE and TM modes in a rectangular waveguide 
filled with isotropic media may be obtained by setting: 


> _ n [Exkz ~ | vuk, 
AA 4) 
which, in our general framework, leads to the solution 
kk, wuk, 
= —jkzz ~ y z ~ u x 
Poe 1) —weky |" “2| kk, 
k kyk 
WEKy yz (55) 
kk, wuk, 
3 kk —wuk 
jkzz ~ y z ~ u x 
+e C3 weky + 4 ~k,k, 
—wek, —kyk, 


Usage of (55) in (26)-(29) gives the TE’/TM? (transverse to 
the z direction) modal decomposition for a waveguide filled 
with isotropic material [12]. We note that if we instead used 
(53), we would obtain the TE*/TE” alternative mode sets 
described by Harrington [12]. We reiterate that alternative 
mode sets do not exist for the waveguide filled with 
anisotropic magneto-dielectric. 


V. CONCLUSIONS 


In this report, a theoretical framework was developed to 
analyze the rectangular waveguide filled homogeneously with 
un-rotated biaxial magneto-dielectric material. It was found 
that the waveguide supported two mode types with reciprocal 
propagation, which for most of the higher order modes, can be 
classified as hybrid, as they are neither TM nor TE to any of 
the Cartesian coordinate axes. The existence of complex 
conjugate pseudo-evanescent mode pairs, and the non- 
existence of alternative mode sets, clearly distinguishes the 
biaxial magneto-dielectric filling from simpler isotropic 
filling. The results obtained in this paper extend the range of 
validity of the results given by Meng et al., providing the 
complete mode spectrum [7]. We note that for more 
complicated anisotropies, including a magnetic material such 
as that considered in this paper, but misaligned with the 
coordinate axes, that an analytic solution cannot be obtained 
by the same methods. For more complex anisotropies, 
numerical techniques are required, and the results will take the 
form of truncated infinite series to approximate the fields and 
propagation constant of a single mode [13], as opposed to the 
closed form expressions (44-50). We note that our analytical 
solution could be utilized to build confidence in these 
numerical schemes, particularly in their ability to predict the 


unusual features discussed in this paper that often have no 
counterpart in problems involving isotropic media, such as the 
hybrid nature of the higher order modes and the existence of 
complex conjugate mode pairs. The analytical solution also 
provides useful theoretical input into development of 
Nicolson-Ross type measurement techniques for 
characterization of anisotropic behavior exhibited by certain 
microwave magnetic materials such as multi-phase magnetic 
composites, single-crystal ceramic ferrites, and thin-film 
ferromagnetic composites [1-6, 8]. 
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